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ABSTRACT 

We examine inhomogeneous primordial nucleosynthesis for arbitrary distributions / of 
the baryon-to-photon ratio t], in the limit where neither particle diffusion nor gravitational 
collapse is important. By discretizing f{rj) and using linear programming, we show that for 
a set of m constraints on the primordial element abundances, the maximum and minimum 
possible values of f] (the final mean value of t]) are given when f{T]) is a sum of at most 
m + 1 delta functions. Our linear programming results indicate that when / is taken to be 
an arbitrary function, there is no lower bound on fj, while the upper bound is essentially the 
homogeneous upper bound. 



1 



1. INTRODUCTION 



One of the most widely-investigated variations in primordial nucleosynthesis is the pos- 
sibility that the baryon density is inhomogeneous at the time of nucleosynthesis. In the 
simplest such inhomogeneous models, nucleosynthesis is assumed to take place indepen- 
dently in regions of different density, with the final element abundances derived by simply 
mixing the matter from these different regions (Epstein & Petrosian 1975; Barrow & Morgan 
1983; Yang et al. 1984; Copi, Olive. & Schramm 1995, hereafter COS). Such a treatment is 
valid if the dominant fluctuations are of sufficiently low amplitude and on a sufficiently large 
length scale that neither gravitational collapse nor particle diffusion is important. If the 
length-scale of the fluctuations is sufficiently short, then particle diffusion before and dur- 
ing nucleosynthesis can significantly alter the element yields (Applegate, Hogan, & Scherrer 
1987; Alcock, Fuller, & Mathews 1987; for more recent studies, see Thomas, et al. 1994; 
Jedamzik, Fuller, & Mathews 1994). On the other hand, if the fluctuation amplitudes are 
high enough, the highest density regions can undergo gravitational collapse, effectively re- 
moving the elements in those regions from the final measured abundances. This scenario has 
also been extensively explored (Sale & Mathews 1986; Gnedin, Ostriker, & Rees 1995; COS; 
Jedamzik & Fuller 1995). 

Here we consider only the simplest of these scenarios (no diffusion or collapse), but we 
provide a significant generalization to previous work. All previous studies have assumed a 
particular form for the distribution of baryon-to-photon ratio f{f])- In this paper, we treat 
the case of arbitrary f{ri). Obviously, a brute-force calculation of this case, examining all 
possible functions, would be impossible. Even the analysis of a limited subset of the infinite 
number of possibilities is an enormous computational task. However, by discretizing the 
problem, we can turn it into a problem in linear programming with a well-defined set of 
easily-obtained solutions. In the next section, we explain our linear programming method 
and calculate the density distributions which give the minimum and maximum possible values 
for the mean value of 77 for a reasonable set of observational constraints on the primordial 
element abundances. Our conclusions are discussed in §3. We find that there is no lower 
bound on 77 (consistent with the earlier arguments of Jedamzik, Mathews, & Fuller 1995), 
while the upper bound is essentially the homogeneous upper bound. 

2. CALCULATIONS 

The most recent and complete treatment of inhomogeneous nucleosynthesis when post- 
nucleosynthesis mixing is the only important effect has been given by COS, so we follow their 
treatment. The elements produced in detectable quantities in primordial nucleosynthesis are 
^He, D, ^He, and ^Li. The observational limits on the abundances of these elements have 
been discussed in great detail in a number of recent papers (Walker et al. 1991; Copi, 
Schramm, & Turner 1995a,b; Hata, et al. 1995, 1996). However, since we wish to compare 
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our results to the calculations of COS, we will use the abundance limits given there: 

0.221 < Yp < 0.243, 
(D/H) > 1.8 X 10-^ 
(D+3He)/H < 1.0 X 10"^ (1) 

COS used two different bounds on the ^Li abundance: 

(^Li/H) < 1.4 X 10"^°, 

(^Li/H) < 2.0 X 10"^°, (2) 

Although we will use these limits in all of our calculations, our method of calculation, which 
is the main new idea in this paper, can be applied to any set of limits. 

Now assume that the distribution of densities is given by an unknown function /(//). 
Previous studies have assumed a variety of functions for f{r}), including the gamma distribu- 
tion (Epstein & Petrosian 1975; Yang et al. 1984) and the lognormal distribution (Barrow 
& Morgan 1983). COS considered both of these as well as a Gaussian distribution for rj. 
We consider here the case of arbitrary /(?]). The mean final value for the baryon-to-photon 
ratio, f], is given by 

roo 

V = / fiv) V dr] (3) 







The final element abundances are mass-weighted averages of the element abundances pro- 
duced in the individual regions. If Xa is the mean mass fraction of nuclide A measured 
today, then 

_ roo 

Xa = / Xa{v) fiv) V dv/v, (4) 







where XA{ri) is the mass fraction of nuclide A produced in standard (homogeneous) nucle- 
osynthesis with baryon-to-photon ratio rj. Since / represents a probability distribution, it 
must be normalized: 

/■oo 

fiv) dv = 1. (5) 







Suppose that we wanted to test all possible functions /. One way to do this would be 
to divide the range in rj into discrete bins, and to approximate the integrals in equations (0) 
and (^) as sums: 

V = J2 fiVi^Vi (6) 

i 

and 

XA = Y.XA^f^V^^Vi/V (7) 

i 

where the rj dependence of / and Xa is expressed through their dependence on the bin 
number i. The bins need not all be of equal size, hence the factor A?7j. The normalization 
constraint (equation |^) becomes 

EM^^ = 1- (8) 
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Now one could imagine doing a Monte Carlo simulation, scanning through all possible 
distributions fi which satisfy the constraint in equation (H). In practice, this is impossi- 
ble. For example, if we divided the range in rj into 100 bins, and divided fi into 1000 
"units" of magnitude 0.001 to be distributed among these bins, we would have to calculate 
C(1099, 1000) ~ 10^^^ sets of element abundances, a calculation clearly beyond all but the 
most robust graduate students. However, most discussions of inhomogeneous nucleosynthesis 
center on a much simpler problem: given a set of constraints such as those given in equations 
(|l|)-(0), and a function or family of functions / for the distribution of f], what are the largest 
and smallest allowed values for f/? If we express the problem this way, discretizing it as in 
equations (|^) and (0), then the question of maximizing or minimizing the value of r] is re- 
duced to a problem in linear programming. (A discussion of both the theory and practice of 
linear programming can be found in, for example. Press, Flannery, Teukolsky, & Vetterling 
1992, from which some of the following discussion is taken). 

The fundamental problem in linear programming is the following: given a set of non- 
negative independent variables Xj, and a set of M constraints of the form: 

N 

^ ^' (9) 

N 

^j^j ^ (10) 

N 

Y,ajXj = b, (11) 

N 

z = J2 CjXj (12) 



or 



or 



maximize or minimize the function 



The linear programming nature of our discretized problem above becomes clearer if we define 
the quantities pt to be 

Pi = fiAr]i. (13) 

Then equations (§)-(!) become 

V = YPiVi, (14) 

i 

XA = Y^AiPiVi/v, (15) 



and 



Ep^ = 1- (16) 



Equations (|14D - (|T6D are now clearly in the form of a linear programming problem: the N 
independent variables are the pj's, equations (|T5|) and (|1^) provide the constraint equations. 
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and fj given by equation (|1^) is the quantity we seek to maximize or minimize. Note that 
equation (|T5|) contains f/ in the denominator, so that the bounds Xa < -^upperbound and 
Xa > -^lowerbound do uot immediately translate into linear programming constraints of the 
form given by equations (pD-(p!T]). However, using equations ([T^ ) and (p!5|), we can rewrite 
these bounds as 

^[^Ai — ^upperbound]Pi^i < 0, (17) 
i 

and 

y^J^Ai — -^lowerbound]Pi^i > 0, (18) 



which are in the form of equations (^ and ([10|) . 



An additional complication is the fact that our constraints on all of the elements other 
than '^He are expressed in terms of number ratios to hydrogen, (A/H)= ua/uh, rather than 
as mass fractions, Xa, while our prescription for mixing the various element abundances uses 
the mass fractions. However, it is easy to translate the number ratio bounds into a suitable 
form. Recall that 

(A/H) = ^, (19) 

where Xh is the ^H mass fraction. Then an observational upper bound on (A/H) of the 
form (A/H) < (A/H) upperbound can be written in the form 

-^A 

^ (^-^ / H) upperbound (^0) 



AX 



H 



where both Xa and Xh are given by equation (pIBf) . Substituting for Xa and Xh from 
equation ([15|) , we obtain 

J2[XAi - (A/H)„boundAXj^i]p,r7, < 0. (21) 

i 

Note that the expression for f/ has dropped out of the equation, but now the inequality 
includes a sum over Xni- For the lower bound (A/H) > (A/H)iowerbound we obtain a similar 
result: 

Yl^Ai - (A/H)iowcrboundAXHi]pir/i > 0. (22) 



Equations (^) and (0) are both in the form of acceptable linear programming constraints. 

We now note an important result of linear programming theory: given a set of X variables 
and M constraints, the solution which maximizes or minimizes z has at least — M of 
the variables Xj equal to zero (Press et al. 1992). In our case, if we have m constraints 
on the element abundances (e.g., equations (|l|) and (^ give m = 5 constraints), and the 
normalization condition (equation ^ provides one additional constraint, then at most M = 
m + 1 of the Pi's are non-zero; in this case M = 6. If we take the continuum limit of equations 
(O)-(H), we arrive at the central result of this paper: for a set of m constraints on the element 
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abundances such as those given in equations (|1|) and (g), and an arbitrary distribution of 77 
given by the function /(r/), the largest and smallest possible values for fj occur when f{ri) is 
the sum of at most m + 1 delta functions. 

We use the simplex method (Press et al. 1992) to solve our discretized problem. We 
ran the primordial nucleosynthesis code of Wagoner (Wagoner, Fowler, & Hoyle 1967) as 
updated by Kawano (1992), with a neutron lifetime of r = 887 sec. [This differs slightly 
from the neutron lifetime in COS, and we do not include the small correction factor to the 
^He production used by COS, but these are small differences which will not significantly 
affect our results]. We allowed 77 to vary from 10~^^ to 10"'', and we divided this interval 
into 600 equally-spaced logarithmic bins, calculating the mass fractions of '^He, ^He, D, and 
''Li for all 600 values of rj. We also included a bin corresponding to zero baryon density 
and zero element production. We then use the simplex method to determine the form for pi 
which maximimizes or minimizes f] for the abundance constraints given in equations (|l|) and 

(!)■ 

When we attempt to minimize t], we find that almost all of pi is concentrated in the 
lowest bin (i.e., the bin corresponding to zero baryon density). In fact, there is no "lowest 
value" for fj; the mean baryon-to-photon ratio can be arbitrarily small, a point recently 
emphasized by Jedamzik, Mathews, and Fuller (1995). For example, one could take /(//) = 
Pi6{ri — rjo) +p2S{ri), where rjo is a value for rj which gives acceptable nucleosynthesis for the 
homogeneous case. By mixing the correct homogeneous rj with the baryon-free regions, we 
obtain the correct element abundances regardless of the values of pi and p2, but by taking 
the limit pi — ^ 0, p2 — ^ 1, the value for f] can be made arbitrarily small. 

A more interesting question is the upper bound on f/. Using either lithium bound in 
equation (0), we find only two non-zero bins, straddling the upper bound on rj in the homo- 
geneous model. To resolve this function further, we re-ran the code using 1000 bins between 
1] = 10~^° and 1] = 10~^. Using the first lithium bound in equation (§) we obtain: 

f{7]) = {0.21)6{r] - 3.31 x 10"^°) + {0.79)6{r] - 3.33 x 10"^°), (23) 

which gives f] = 3.33 x 10~^°. In this case, it is the limit on ^Li which is saturated. The 
second lithium bound in equation (0) yields the solution: 

f{r]) = (0.58)5(7] - 3.40 x 10"^°) + (0.42)5(r/ - 3.44 x 10"^°), (24) 

with fj = 3.42 X lO"^''. For this case, the ^He limit is saturated. In fact, given the precision 
with which we calculate the various element abundances, the difference between these bounds 
and the homogeneous upper bound is not significant. We are led to conclude that the 
homogeneous upper bound on r] cannot be exceeded for any distribution /(r/). Our results 
are consistent with the claim that /(t/) should be at most a sum of m + 1 delta functions; in 
this case m + 1 = 6, while our optimum solution is the sum of only two delta functions. 

Another interesting question is whether the inclusion of inhomogeneities can resolve a 
slight discepancy between the predictions of primordial nucleosynthesis and the observations. 
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A number of recent studies (Copi, Schram, & Turner 1995a,b; Hata, et al. 1995) have 
suggested that values of r) which give production of D, ^He, and ^Li consistent with the 
observations will overproduce ^He, so it has been suggested, for example, that the "true" 
primordial ^He abundance is somewhat larger than is currently reported by the observers. 
Since inhomogeneities allow for the reduction of fj by an arbitrary amount, is it tempting to 
think that a reduction in the ^He abundance might also be possible. To test this, we have 
used the D, ^He, and ^Li limits given above, while minimizing the value of ^He. Our linear 
programming method cannot be applied directly to minimize ^He, because the expression for 
^He contains a factor proportional to fj. However, we can find the distribution in pi which 
minimizes ^He for any given value of f], so we have simply scanned over a range in fj to find 
the smallest ^He as a function of fj. Using this method, we find that no significant reduction 
in ^He is possible. Again, our results apply to arbitrary distributions of rj. 

3. DISCUSSION 

We find that for this simplest model for inhomogeneous nucleosynthesis, there is no lower 
bound on 77, as expected (Jedamzik, Mathews, & Puller 1995). Previous studies which as- 
sumed particular functional forms for f{r)) all produced a fairly narrow range in the allowed 
values for fj (Epstein & Petrosian 1975; Barrow & Morgan 1983; Yang et al. 1984; COS). 
However, this occurred because all of the functions /(?]) in these papers were unimodal, i.e., 
characterized by a single maximum. Hence, they cannot approximate the sort of solution 
which has two large peaks at ?] = and rj = rjo (where rjo gives acceptable abundances in the 
homogeneous model). Our upper bound on fj is, for all practical purposes, the homogeneous 
upper bound on rj. This is consistent with the results of COS; for all three functions they 
examined, it is clear from their figures that the upper bound on fj is no larger than that 
obtained when the variance of /(?]) goes to zero. The importance of our results is that they 
give the most general upper and lower bounds on fj for any density distribution, essentially 
bringing to a close the two-decade-long investigations of these simplest inhomogeneous mod- 
els. Our final conclusion is that for the inhomogeneous models in which mixing is the only 
important process, fj can be arbitrarily small (a point already noted by Jedamzik, Mathews, 
& Puller 1995), but fj cannot be larger than the homogeneous upper bound on rj. 

Linear programming is not a technique usually applied to astrophysical problems, al- 
though it has been previously used in galactic dynamics (Schwarzschild 1979). Our technique 
of discretizing the problem and using linear programming could be applied to any problem 
with constraints on integrals of an unknown function. This technique cannot be applied to 
inhomogeneous nucleosynthesis when particle diffusion is significant. However, it could be 
applied to the case when collapse of high density regions is important. We have chosen not to 
address this case because it involves many particular model assumptions, but, for example, 
equation (21) of COS can easily be put into the form of a linear programming constraint. 
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